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Abstract 

We survey the non-locked color-flavor-spin channels for quark-quark (color 
superconducting) condensates in QCD, using an NJL model. We also 
study isotropic quark-antiquark (mesonic) condensates. We make mean-field 
estimates of the strength and sign of the self-interaction of each condensate, 
using four-fermion interaction vertices based on known QCD interactions. For 
the attractive quark pairing channels, we solve the mean-field gap equations 
to obtain the size of the gap as a function of quark density. We also calculate 
the dispersion relations for the quasiquarks, in order to see how fully gapped 
the spectrum of fermionic excitations will be. We use our results to specify 
the likely pairing patterns in neutral quark matter, and comment on possible 
phenomenological consequences. 



1 Introduction 



It is by now well known that the BCS mechanism that underlies superconductivity 
in metals is likely to operate even more strongly in dense quark matter [1, 2, 3, 4, 5] 
(for reviews, see Ref. [6]). Generally, pairing between quarks of different flavors 
has received most attention. This is because the strongest attractive interaction 
occurs in the color-antisymmetric spin-zero channel, so by Fermi statistics the flavor 
wavefunction must be antisymmetric, involving two different flavors. This is the 
pattern of pairing in the heavily-studied "2SC" [4, 5] and "CFL" [7] phases of quark 
matter. However, there is strong reason to believe that these are not the only pairing 
patterns that are relevant in nature, and in this paper we survey and discuss some 
of the channels that are available to quarks that cannot participate in 2SC or CFL 
pairing. 

The essence of the BCS mechanism is that a Fermi surface is unstable against 
pairing if there are any fermion-fermion channels in which the interaction is 
attractive [8]. For electrons in metals this condition is only met where phonon- 
mediated attraction overwhelms Coulomb repulsion. But there are certainly 
channels in which the dominant QCD quark-quark interaction is attractive, as 
indicated by the fact that at low density and temperature quarks bind strongly 
together to form hadrons. Mean-field calculations using QCD-inspired NJL models 
[7, 9, 10, 11], lattice studies of NJL models [12], and calculations using gluon 
exchange with a hard dense loop resummed gluon propagator [13, 14, 15, 16, 17] 
confirm this, and indicate that the main pairing patterns are CFL (for sufficiently 
light strange quark or high density), or 2SC/2SC+S for heavier strange quarks. 
However, these calculations do not take into account the requirement of electric and 
color neutrality, which must be obeyed in uniform dense matter in the real world, 
such as might be found in the core of a compact star. 

The effects of imposing neutrality were recently studied in a model-independent 
expansion in powers of m s /fi [18], and in preliminary Nambu-Jona-Lasinio model 
calculations [19]. Those calculations indicate that the 2SC+S phase suffers a 
significant penalty for remaining neutral, and may not occupy such a large part 
of the phase diagram as previously thought. 

A conjectured phase diagram is sketched in Fig. 1, where we show schematically 
the regions of phase space in which various phases, each breaking a particular 
set of symmetries, are expected to exist. The diagram shows the m s -fi plane at 
temperature T = 0, and the u and d quark masses are set to zero. 

At low strange quark mass (m s <C m c ° nt ), compression of nuclear matter leads 
to the production of hyperons (the "strange hadronic phase"), and then a transition 
into an isospin-symmetric phase which by quark-hadron continuity [20] can be 
interpreted as a particular pairing pattern of baryons or as quark matter in the 
chiral symmetry breaking color-flavor-locked phase. 

At slightly higher strange quark mass, compression of strange hadronic matter 
may lead into an isospin broken phase that can be interpreted as hyperonic matter 
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Figure 1: A conjectured phase diagram for neutral dense matter at zero temperature 
as a function of quark chemical potential \i and strange quark mass m s according 
to Ref. [18]. The up and down quark masses are assumed zero. The symmetries 
broken by the various phases are indicated by the different shading styles. The 
phase structure in the chirally unbroken wedge in the upper right hand corner is 
still uncertain, and is the topic of this paper. 

with S° pairing [21] or as color-flavor-locked quark matter with K° condensation 
( "CFL-K " ) [22] . However, the existence of this phase at these densities is uncertain 
because of possibly large instanton effects [23] . 

At high strange quark mass (m s > m c s ont ) the CFL phase is pushed to higher 
densities, and there is an interval of densities where we expect non-color-flavor-locked 
quark matter to exist. This is the wedge in the upper right part of the m s -fi plane 
(Fig. 1) where chiral symmetry is restored. (In the figure it is assumed that close to 
the transition to CFL there will be some strange quarks present, and they will pair 
somehow, breaking strangeness.) We are uncertain of the favored phase(s) in that 
region. Since only the CFL phase pairs all the quarks that are present, this region 
will involve phases that leave single "orphaned" flavors or colors of quark to find an 
attractive channel in which to pair. We can distinguish two general scenarios: 
(1) There may be no BCS pairing of different flavors [18], so all nine quark 
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color/flavor species are left to find alternative pairing channels. 
(2) As the calculation of Ref. [19] suggests, there may be some 2SC+S, i.e. a phase 
that pairs the u and d quarks of two colors (red and green, say), leaving all the s 
quarks and the blue u and d quarks unpaired. 

In both cases the BCS argument implies that the supposedly unpaired quarks will 
in reality seek some attractive channel in which to pair, and in this paper we will 
study some of the possibilities. This question is of direct relevance to compact star 
physics, where the smallest gap controls transport properties such as the specific 
heat and neutrino emission rate. The density of a compact star rises from nuclear 
density near the surface to a much higher but as yet unknown value in the core, so 
unless there is a direct transition from nuclear to CFL quark matter, there will be 
regions where some sort of single color or flavor pairing is likely. 

One possibility for the orphaned quarks is crystalline or "LOFF" pairing [24, 25, 
26] in which species whose Fermi surfaces are too far apart to support standard 
translationally-invariant BCS pairing instead form pairs with net momentum, 
utilizing only part of their Fermi surfaces, and breaking translational invariance. In 
this paper we will be concerned with another possibility: single-color or single-flavor 
translationally invariant (but not necessarily isotropic) pairing [3, 4, 27, 28]. Such 
phases may have non-zero angular momentum, spontaneously breaking rotational 
invariance. In the following sections we survey the non-locked quark pairing 
channels, paying particular attention to the single flavor and/or single color channels. 
We identify the attractive channels using NJL models with four-fermion interactions 
based on instantons, magnetic gluons, and combined electric and magnetic gluons. 
For the attractive channels we solve the gap equations, deduce the quasiquark 
spectrum, and comment on possible physical manifestations. Finally, we perform a 
similar survey of mesonic channels, in which exotic pairing has been posited [29]. 



2 Mean-field survey of quark pairing channels 
2.1 Calculation 

To see which channels are attractive we perform a mean-field calculation of the 
pairing energy for a wide range of condensation patterns. We write the NJL 
Hamiltonian in the form 

H Hf Tee -\- H interaction 

H {iee = - fir/o + m)i/j (2.1) 

j ^fl n-ia f3b kc 5d > 



interaction = KVf^Vf W" 



where color indices are a, j3, 7, 5, flavor indices are i, j, k, I, spinor indices are a, b, c, d. 
The four-fermion interaction is supposed to be a plausible model of QCD, so in the 
interaction kernel H, we include three terms, with the color-flavor-spinor structure 
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of a two-flavor instanton, electric gluon exchange, and magnetic gluon exchange, 

^clcc 7~~tmag ^inst 

H elec = IG E 5i5[ 5 ab 5 cd 1(35^5} - 6%6]) 

H mag = \G M 5\5{ ELl[707n]a b [707n]e, §(3^ - 5%6]) 

^inst = e^e*' |([7o(l + 7s)]afe[7o(l + lb)]cd + [7o(l - 7b)U[7o(1 - lh)]cd) 

1(35^5] -5f5j) 

(2.2) 

We consider condensates that factorize into separate color, flavor, and Dirac 
tensors (i.e. that do not show "locking") and calculate their binding energy by 
contracting them with (2.2). 

There is no Fierzing ambiguity in this procedure. For a given pairing pattern 
X, the condensate is 

(VffyfV, = A(X)€^ (x)ji r b { d x) . (2.3) 
We can then calculate the interaction ( "binding" ) energy of the various condensates, 

# = &( x ) 2 ( s a2 G E + SSHGm + sffiG,) (2.4) 

fx) 

The binding strengths interaction & ve the strength of the self-interaction of the 
condensate X due to the specified part of the interaction Hamiltonian. 

2.2 Properties of the pairing channels 

In Table 1 we list the the simple (translationally invariant, factorizable) channels 
available for quark pairing. The meanings of the columns are as follows. 

1. Color : two quarks either make an antisymmetric color triplet (which requires 
quarks of two different colors) or a symmetric sextet (which can occur with 
quarks of two different colors, and also if both quarks have the same color). 
For the 3a we use £P S = e 136 in Eq. (2.3). For the 6s we use a single-color 
representative & & = <5 /3 ' 1 <5 5 ' 1 in Eq. (2.3). 

2. Flavor : two quarks either make an antisymmetric flavor singlet (which requires 
quarks of two different flavors) or a symmetric triplet (which can occur with 
quarks of two different flavors, and also if both quarks have the same flavor). 
For the 1a we use 3^ = &\ and for the 3s we use 3.?; = aj t in (2.3). 

3. Spin, parity: since the chemical potential explicitly breaks the Lorentz group 
down to three-dimensional rotations and translations, it makes sense to classify 
condensates by their total angular momentum quantum number j and parity. 
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4. Dirac : This column gives the Dirac matrix structure T bd used in (2.3), so the 
condensate is ip T Ti(j. We also designate each condensate as "LL" (even number 
of gamma matrices, so pairs same-chirality quarks) or "LR" (odd number of 
gamma matrices, so pairs opposite-chirality quarks). 

5. BCS-enhancement : Condensates that correspond to pairs of particles or holes 
near the Fermi surface have a BCS singularity in their gap equation that 
guarantees a solution, no matter how weak the coupling. To see which 
condensates have such a BCS enhancement, we expanded the field operators 
in terms of creation and annihilation operators (see Appendix C). The order 
of the coefficient of the a(p)a(— p) and tf(p)tf(— p) terms is given in the 
table. 0(1) means BCS-enhanced, means not BCS-enhanced. In the C7075 
condensate the coefficient goes to zero as the quark mass goes to zero (hence 
it is labelled "(9(m)" in the table) meaning that the channel loses its BCS 
enhancement in the chiral limit. This is discussed further in Appendix C. 

6. Binding strength: For each channel we show the binding strength for the 
instanton interaction, the full (electric plus magnetic) gluon, which could 
reasonably be used at medium density, and for the magnetic gluon alone, which 
is known to dominate at ultra-high density [30, 13]. Channels with a positive 
binding strength and BCS enhancement will always support pairing (the gap 
equation always has a solution, however weak the couplings Gi,Ge,Gm)- 
Other things being equal, the pairing with the largest binding strength will 
have the lowest free energy, and is the one that will actually occur. 

It may seem strange that there are entries in the table with angular momentum 
j = 1 and an antisymmetric Dirac structure (C7375), and with j = but a 
symmetric Dirac structure (C70). If all the angular momentum came from spin this 
would be impossible. But even though there are no explicit spatial derivatives in the 
diquark operators, there can still be orbital angular momentum. In Appendix C the 
angular momentum content of the particle-particle component of the condensates 
is analyzed into its spin and orbital content. We see, for example, that C7375 has 
an antisymmetric space wavefunction (I = 1) and a symmetric spin wavefunction 
(s = 1), combined to give an antisymmetric j = 1. 

2.3 Results 

The results of the binding strength calculation are shown in Table 1. The first block 
is antisymmetric in flavor and color, and so describes pairing of two flavors and two 
colors. The second block is for two flavors and one color, the third for one flavor 
and two colors, and the final block is for one color and one flavor. 

Certain features can be easily understood: the flavor-symmetric condensates 
all have zero instanton binding energy, because the instanton vertex is flavor- 
antisymmetric in the incoming quarks. The gluonic vertices give the same results for 
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Table 1: Binding strengths of diquark channels in NJL models in the mean-field 
approximation. The first 6 columns specify the channels, and the last 3 columns 
give their attractiveness in NJL models with various types of four-fermion vertex: 
2-flavor instanton, single gluon exchange, single magnetic gluon exchange (expected 
to dominate at higher density). See equations (2.3) and (2.4) and subsequent 
explanation. 

C75 as for C, and for C 00375 as for CV03, because the gluonic interaction is invariant 
under U(1)a transformations, under which the LL condensates transform into each 
other (C75 C and Ccr 375 ^ Ca 03 ) while the LR condensates are invariant. We 
see that there are many attractive channels: 
1) Two colors and two flavors (3a, 1a,- • •)• 

The strongly attractive channel (3a,1a,0,+)(C7s) is the 2SC and CFL quark 
Cooper pairing pattern, and has been extensively studied. The gap is large 
enough that even species with different masses, whose Fermi momenta are 
quite far apart, can pair (hence the CFL phase which pairs red and green u and 
d, red and blue u and s, and green and blue d and s in this channel). Its parity 
partner (3a,1a,0,— )(C) is disfavored by instantons, and is therefore unlikely 
to occur at phenomenologically interesting densities. The additional channel 
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(3^,1^,1,— )(C*7375) is more weakly attractive and also breaks rotational 
invariance, and is therefore expected to be even less favored. This is confirmed 
by gap equation calculations (Fig. 3) which show that its gap is smaller by a 
factor of 10 to 100. 

2) One color, two flavors (6s, 1a,- ■ •)• 

It is generally emphasized that the quark-quark interaction is attractive in 
the co lor- antisymmetric 3 a channel. But, as we see in table 1, the color- 
symmetric (65,l J 4,l,+)(Co"o3) is attractive for instantons and the magnetic 
gluon four-fermion interaction. The instanton gives it a gap of order 1 MeV 
(Fig. 2), while the gluon interaction gives a small gap of order 1 eV (Fig. 3). 
This channel was originally suggested for pairing of the blue up and down 
quarks that are left out of 2SC [4], and is discussed in more detail in Ref. [28]. 
Its gap is small, so it could only pair quarks of similar mass, i.e. the light 
quarks, but in a real-world uniform phase such pairing will not occur either, 
because charge neutrality causes the up and down chemical potentials (and 
hence Fermi momenta) to differ by tens of MeV, which is larger than the gap. 
In a non-uniform mixture of two locally charged phases [31], however, it is 
conceivable that the up and down Fermi momenta could be similar enough 
to allow pairing in this channel. The parity partner (65,1^,1,— )(Ca^ r )^) is 
disfavored by instantons. The channel (6s,1a,0,— )(C*7o) is attractive, but has 
no particle-particle component, and presumably only occurs for sufficiently 
strong coupling. Solving the gap equations for reasonable coupling strength 
we find no gap in this channel. 

3) Two colors and one flavor (3^, 3s,. . .). 

The only attractive channel is (3^,35,1,+) (C73). This is a pairing option for 
red and green strange quarks in 2SC+S. We have solved the relevant gap 
equation (Figs. 3,4) and find gaps in the range 2 — 10 MeV. If three colors 
are available then a competing possibility is to lock the colors to the spin 
(CSL), so the condensate is a linear combination of and Ca^ with a 
color structure that is correlated with the spatial direction, e.g. red and green 
quarks pair in the z direction, red and blue in the y direction, green and blue 
in the x direction. This leaves an unbroken global 50(3) of spatial rotations 
combined with color rotations, so the gap is isotropic, which helps to lower 
the free energy [27, 32]. Note also that the channels (3a,3s,1,+)(Cct 03 ) and 
(3^,35,1,— )(Ccr 375) which are repulsive in the NJL model become attractive 
at asymptotic density when the gluon propagator provides a form factor that 
strongly emphasizes small- angle scattering [27]. 

4) One color and one flavor (65,35-,. . .). 

There is an attractive channel here, the (65,35,0,+) (C7075). It loses its 
particle-particle component as the quark mass goes to zero, making it very 
weak for up and down quarks, but quite strong for strange quarks (Fig. 4). It 
is suitable for the blue strange quarks in 2SC+S when red and green strange 
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quarks have paired in the (3^,3,5, 1,+)(C7 3 ) channel. 
Many of the attractive channels have repulsive partners with the same 
symmetries, so a condensate in the attractive channel will automatically generate a 
small additional one in the repulsive channel. For example, the (3^,1,4,0,+) (C75) 
can generate (3 J 4,l J 4,0,+)(C , 7o75). This was discussed in Ref. [33], where the induced 
(C7075) condensate (there called "k") in 2+1 flavor CFL was calculated and found 
to be small. In Ref. [34] it was observed that if all three quarks are massive then 
this condensate may be important. In the context of CFL the (3a,1a,0,+)(C75) 
can also generate (6s,3s,0,+)(C'j5) [7, 35], since they both break the full symmetry 
group down to the same subgroup. 

Instanton Interaction 
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Figure 2: Gap parameters in the attractive channels for an NJL interaction based on 
the two-flavor instanton. Since the instanton interaction requires two quark flavors, 
we take the quarks to be massless, which is a good approximation for the u and d. 
The cutoff is A = 800 MeV. 



3 Gap calculations for the attractive diquark channels 

For the attractive channels we performed uncoupled gap equation calculations, and 
obtained the dependence of the quark pairing on \x. The amount of pairing is given 
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Figure 3: Gap parameters in the attractive channels for an NJL interaction based 
on magnetic-gluon exchange. We show the one-flavor and two-flavor channels, for 
massless quarks. The cutoff is A = 800 MeV. 

by the gap parameter A(/i), which occurs in the self energy (See Appendix A) as 

A°f a6 (p) = A(/x) C a ^ ij T ab , (3.1) 

with color matrix C, flavor matrix J 7 , and Dirac structure T ab . Note that A(/i) 
is a gap parameter, not the gap. It sets the scale of the gap in the quasiparticle 
excitation spectrum, but as we will see in Section 4 the gap itself often depends on 
the direction of the momentum. 

The 4-fermion interactions that we use are nonrenormalizable, so our gap 
equation involves a 3-momentum cutoff A, which represents the decoupling of our 
interactions at higher momentum, due to instanton form factors, effective gluon 
masses, etc. The usual procedure for NJL model calculations is to calibrate the 
coupling strength for each cutoff A by known low-density physics such as the size of 
the chiral condensate. However, it is well known that this leads to an approximately 
cutoff-independent maximum gap (as a function of n) in the ipC^ijj channel, so we 
used that criterion directly as our calibration condition, setting the maximum gap 
to 100 MeV. 
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Figure 4: Gap parameters in the attractive channels for an NJL interaction based on 
magnetic-gluon exchange, for quarks of mass 250 and 350 MeV, a reasonable range 
of values for the strange quark at medium density. We show only the single-flavor 
channels. The cutoff is A = 800 MeV. 

The results of our calculations, for cutoff A = 800 MeV, are plotted in Figs. 2, 
3,4. For other cutoffs the overall shape of the curves is very similar. Because we use 
a sharp cutoff A, the gap falls to zero when \i reaches A (see, e.g., Eq. (A. 16)). We 
show gap plots for the instanton interaction (Fig. 2) and magnetic gluon interaction. 
The full electric + magnetic gluon gives results that are similar to those for the 
magnetic gluon, but with no gap in the (65,1,4,1,+) (Ca 03 ) channel. 

For the magnetic gluon, we show a gap plot for massless quarks (Fig. 3) 
which includes the two-flavor channels (3a,1a,0,+)(C7s), (3^,1^,1,— XC7375) and 
(65,1^,1,+) (C0O3) which could sustain u-d pairing, as well as the single-flavor 
channel (3a,3s,1,+)(C73) which could sustain u-u or d-d pairing. We also show 
a gap plot for quarks with a mass appropriate to the strange quark (Fig. 4) which 
includes the single-flavor channels that could sustain s-s pairing with two colors 
(3 A ,3 5 ,1,+)(C , 7 3 ) or one color (6 S ,3 5 ,0,+)(C7 7 5 ). 

The relative sizes of the gaps in the different channels reflect the pairing strengths 
given in Table 1. We see that the Lorentz scalar (3^,1^,0,+) (C75) (red/solid 
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line) is dominant. The j — 1 channels have much smaller gap parameters. The 
(3^4,35,1,+) (C73) gap parameter (blue/dash-dot line in Figs. 3 and 4) rises to a few 
MeV with the magnetic gluon interaction. The (6s,l^,l;+)(C cr 03) gap parameter 
(green/dashed line) rises to about 1 MeV with an instanton interaction, but only 
1 eV with the magnetic gluon interaction. It should be remembered, however, that 
the temperature of a compact star can be anything from tens of MeV at the time 
of the supernova to a few eV after millions of years, so gaps anywhere in this range 
are of potential phenomenological interest. 

The (65,35,0,+) (C7075) channel (black/dotted line), which is the only attractive 
channel for a single color and flavor of quark, is highly suppressed for massless quarks 
at high density but reaches about 10 keV for strange quarks (m = 350 MeV, Fig. 4). 
This is because its particle-particle component goes to zero as m — > (Eq. (C.9) 
and Table 1). 

Up to this point we have not mentioned the j = l,mj = ±1 channels 
(e.g. ipC^±ip = ipC^i ±272)^)- We have only discussed the j = l,rrij = channels 
(e.g. i/jC'jsip). That is because rotational invariance of the interaction Hamiltonian 
that we are using guarantees that changing mj from to ±1 will not affect the 
binding energy and gap equation. This can be seen by considering the form of the 
binding energy. From Eq. (2.1) it is 

E B ~ (^) ]ac (^) bd H abcd . (3.2) 

Note that it is quadratic in the diquark condensate, with one of the factors 
being complex conjugated. So if we have some 3-vector condensate, for example 
4> = J2i 0i(V , 7i'0)5 then its binding energy is 

£ B cx|^.| 2 + |</g 2 + |</g 2 (3.3) 

It is clear that the mj = condensate ^ = (0, 0, 1) has the same binding energy 
as the mj = ±1 condensate 0j = (l/v2)(l, ±z, 0). We have explicitly solved the 
gap equations for the mj = ±1 condensates, and find their solutions identical to the 
corresponding mj = condensates. However, the quasiquark excitations in the two 
cases are quite different, and we proceed to study these in the next section. 

4 Quasiquark Dispersion relations 

The physical behavior of quark matter will be dominated by its lowest energy 
excitations. As well as Goldstone bosons that arise from spontaneous breaking of 
global symmetries, there will be fermionic excitations of the quarks around the Fermi 
surface. In the presence of a diquark condensate, the spectrum of quark excitations 
is radically altered. Instead of arbitrarily low energy degrees of freedom, associated 
with the promotion of a quark from a state just below the Fermi surface to just 
above it, there is a minimum excitation energy (gap), above which the excitation 
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spectrum is that of free quasiquarks, which are linear combinations of a particle and 
a hole. 

The dispersion relations of the quasiparticles can be calculated straightforwardly 
by including a condensate of the desired structure in the inverse propagator S*" 1 , 
shown in Eq. (A.3). Poles in the propagator correspond to zeros in S*" 1 , so the 
dispersion relations are obtained by solving det S~ 1 (pq, p, /i, A, m) = for the energy 
Po as a function of the 3-momentum p of the quasiparticle, quark chemical potential 
/i, gap parameter A, and quark mass m. 

The gap is by definition the energy required to excite the lowest energy 
quasiquark mode. Isotropic condensates have a uniform gap, but one of the most 
interesting features of j > condensates is that they are not in general fully gapped: 
the gap goes to zero for particular values of momentum p, which correspond to 
particular places on the Fermi surface. This means that transport properties such 
as viscosities and emissivities, which are suppressed by factors of exp(— A/T) in 
phases with isotropic quark pairing, may not be so strongly suppressed by a j > 
condensate. In Figs. 5 and 6 we show the variation of the gap over the Fermi surface 
by plotting the energy of the lowest excitation as a function of angle, 

£ gap (0)=min 1^,0)1 (4.1) 

p.i 

where Ei(p,6) is the energy of the i th quasiquark excitation with momentum 
(psin^cos^psin^sin^pcos^). For the plots we take fi = 500 MeV and A = 
50 MeV, with quark mass m = (Fig. 5) or m = 250 MeV (Fig. 6). 

1) C73 condensate: j = l,rrij = 0. 

There is one quasiquark excitation with energy less than the gap parameter A. 

/ieffW 2 = ^ 2 + A 2 cos 2 (£) (4.2) 
A cS (6) 2 = A 2 (sin 2 (0)+m 2 / yUc 2 ff cos 2 (0)) 

From the expression for A e ff (6) we see that for massless quarks the gap goes to zero 
for momenta parallel to the z-axis, i.e. at the poles on the Fermi surface (red/solid 
curve in Fig. 5). Massive quarks retain a small gap of order mA//j at the poles 
(red/solid curve in Fig. 6). 

2) C(7i ± «7 2 ) condensate: j = l,rrij = ±1. 

There are two quasiquark excitations with energy less than 2 A, 

E(p) 2 = 2A 2 + m 2 + + p 2 ± (4A 4 + 4/i 2 (p 2 + m 2 ) + 2A 2 p 2 (1 - cos(2fl)) 

V i (4.3) 

±4A 2 /i v /4m 2 + 2p 2 (l + cos(2#)) j 2 

For this condensate the effective gap again goes to zero at the poles, but in this 
case it remains zero even in the presence of a quark mass (magenta/dotted curve in 
Figs. 5,6). 
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Figure 5: Energy gap in units of the gap parameter as a function of polar angle on 
the Fermi surface for rotational symmetry breaking phases with massless quarks, at 
/i = 500 MeV, gap parameter A = 50 MeV. 7+ = 71 + ry 2 , cr+ = 001 + ^02- 



3) Ca 03 condensate: j = l,rrij = 0. 

There is one quasiquark excitation with energy less than A, its dispersion relation 
is [28] 

E( P ) 2 = (Vp 2 + mV/^± /lcff ) 2 + A c 2 ff 
/XeffW 2 = ^ 2 + A 2 sin 2 (fl) (4.4) 
A cff (#) 2 = A 2 (cos 2 (#) W//i 2 ff sin 2 (#)) 

This is related to the dispersion relation for the C73 condensate by 9 — > tt/2 — 0: 
for massless quarks the quasiquarks are gapless around the equator of the Fermi 
sphere (blue/dash-double-dot curve in Fig. 5) and in the presence of a quark mass 
they gain a small gap of order Am//x (blue/dash-double-dot curve in Fig. 6). The 
equator is a larger proportion of the Fermi surface than the poles, so in this case we 



13 




Polar angle 9 

Figure 6: Energy gap in units of the gap parameter as a function of polar angle on 
the Fermi surface for rotational symmetry breaking phases at fi = 500 MeV, with 
gap parameter A = 50 MeV. The quarks have mass m = 250 MeV. 7 + = 71 + 272, 

0"+ = CM + «O"02- 

might expect a greater effect on transport properties. 
4) C(<7oi ± condensate: j = l,rrij = ±1. 

There are two quasiquark excitations with energy less than A. They have rather 
complicated dispersion relations. Going to the massless case, and assuming 
E, (p — fi) <C fi, which will be true for the low-energy quasiquark degrees of freedom 
that we are interested in, we find 

E(p) = (AV + AVos(fl) ± r//v / 2)/(2/i 2 ) , . 

= (8/i 4 (/i - pf + 8A 2 /i 2 (Ai 2 - cos(0)V) + 2A 4 (/i + cos(^)p) 2 ) K ' 

In this case there is a region near the poles, 9 < A//i, where the gap is zero 
(green/dashed curve in Fig. 5). This is because at those angles E{p) has zeros at 
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two values of p close to \i. When 9 m A/fj, those two zeros merge and disappear 
from the real p axis. The presence of a quark mass m > A wipes out this effect, 
but there is still no gap at the poles on the Fermi surface (green/dashed curve in 
Fig. 6). 

In comparing Figs. 5 and 6 it is interesting to note that introducing a mass for 
the quark opens up a gap whenever the gap lines intersect each other at a non- 
zero angle (after one includes the mirror-image negative-energy gap curves for the 
quasiholes). This occurs at zero energy at the poles for C73 and at the equator for 
C<To3. It occurs at non-zero energy for C(ji ± 172). The case of C(<7oi ± icr^) is 
similar, but it is not obvious from the gap plot for reasons described above. 

We see that the j ^ phases show a rich variety of quasiquark dispersion 
relations. For massless quarks they are all gapless in special regions of the Fermi 
surface, and for massive quarks the rrij = ±1 condensates remain gapless for 
momenta parallel to the spin. It follows that for these phases the quasiquark 
excitations will play an important role in transport properties, even when the 
temperature is less than the gap parameter. 

Moreover, different condensates (rrij = ±1 vs. rrij = 0) that because of rotational 
invariance of the Hamiltonian have exactly the same binding energy and gap 
equation, nevertheless have completely different energy gaps over the Fermi surface. 
They will therefore behave quite differently when exposed to nonisotropic external 
influences, such as magnetic fields or neutrino fluxes, and also in their coupling 
to external sources of torque, e.g. via electron-quasiquark scattering. All these 
influences are present in compact stars, and it will be interesting and complicated 
to sort out which is favored under naturally occurring conditions. And it should not 
be forgotten that these conditions vary with the age of the star. 

5 Mesonic condensates 

Berges and Wetterich [29] have invoked ideas of complementarity to suggest that 
the confining QCD vacuum could be understood in terms of non-color-singlet chiral 
condensates. We can study the attractiveness of such channels by mean-field 
methods similar to those of section 2. 

A mesonic condensate X with strength Apn is 

^f) = A {x) ^ {x) %V (x) b a (5.1) 

where is the color-flavor structure of the mesonic condensate. We can calculate 
their binding energy by contracting them with (2.2), and the energy is given by 

( X) 

(2.4), where as before the binding strengths £i nterac tion & ve the strength of the 
self-interaction of the condensate X due to the specified part of the interaction 
Hamiltonian. Unlike the diquark condensates, there are two possible contractions, 
the Hartree contribution and the Fock contribution, since there are two possible ip^ 
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to contract with each ip. Again there is no Fierz ambiguity, and we include both 
Hartree and Fock contributions. 

For consistency we have used the same three interactions, namely instanton, full 
(electric plus magnetic) gluon, and magnetic gluon, as in our treatment of quark 
pairing. However, it should be noted that the main relevance of quark- ant iquark 
pairing is at low density, where there is no reason to expect the magnetic gluon to 
predominate, so in the case of Table 2 the final column does not have any special 
physical relevance. 

The color-flavor structures that we study are 



The (1,1) is a singlet in color and flavor. The (8,8) is a color-flavor locked adjoint 
chiral condensate, of the kind posited in Ref. [29]. We have studied it in the two- 
flavor case, where the two flavors lock to two of the three colors (red and green, 
say). The (8,1) is a flavor-singlet condensate with color structure (1,1,-2), so it 
spontaneously selects a color direction that we have fixed as blue. 

The Dirac structures are two chiral condensates 1 and 75 (scalar and pseudoscalar 
respectively), and two that contain a 70 and therefore correspond to a number 
density of quarks. Nr + Nl for 70, i.e. a renormalization of the chemical potential, 
and N R - N L for 7075. 

The results are given in table 2. We see that the standard chiral condensate 
(1, 1, 1) is overwhelmingly favored over all the others. Its pseudoscalar partner 
(1,1) (75) is disfavored by instantons. 

Among the color-flavor-locked adjoint condensates, the interaction energies are 
much weaker. The Lorentz scalar (8,8) (1) studied in Ref. [29] is disfavored by all 
the interactions, but its pseudoscalar partner (8,8) (75) is less disfavored, thanks 
to a positive contribution from instantons. Interestingly, the only favored adjoint 
condensates are the ones that correspond to spontaneous generation of particle 
number, the (8,8) (70) and its pseudoscalar partner. 

The color-adjoint flavor-singlet condensates show larger interaction energies than 
the color-flavor-locked adjoint condensates, but mostly they are repulsive. The very 
large repulsion for the (8,1) (70) in the electric gluon channel comes from the Hartree 
term, and corresponds to the enormous Coulomb repulsion of a spontaneously 
generated uniform density of blue quarks. This color imbalance does not arise in 
other condensates: in the color-flavor locked (8,8) (70) each flavor favors a different 
color, and in the (8,8) (27075) the left-handed and right-handed quarks have opposite 
color. If instantons predominate then the (8,1) (1) condensate would be the strongest 
competitor to the standard chiral condensate. 



m 



(5.2) 
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Structure of 
condensate 

color-flavor Dirac 



Attractiveness 



mag + elec magnetic 
mstanton , f 

gluon gluon 



(1,1) 


1 


+192 


+192 


+ 144 


(1,1) 


75 


-192 


+ 192 


+144 


(1,1) 


7o 





-96 


-144 


(1,1) 


«7o75 





-96 


-144 


(8,8) 


1 


-6 


-12 


-9 


(8,8) 


75 


+6 


-12 


-9 


(8,8) 


7o 





+6 


+9 


(8,8) 


«7o7s 





+6 


+9 


(8,1) 


1 


+24 


-48 


-36 


(8,1) 


75 


-24 


-48 


-36 


(8,1) 


7o 





-300 


+36 


(8,1) 


^7o75 





-12 


+36 



Table 2: Quark- ant iquark pairing strength in various rotationally symmetric 
channels. 

6 Summary 

We have surveyed the factorizable (non-locked) color-flavor-spin quark pairing 
condensates that occur in an NJL model of QCD that uses point-like four-fermion 
interactions (Table 1). We have solved the gap equation for the attractive channels, 
and obtained the relative sizes of their gap parameters (Figs. 2,3,4). 

We conclude that there are three types of neutral quark matter that may be 
found in the mysterious non-CFL wedge in the high strange quark mass and high 
density region of the m s -[i plane (Fig. 1): two-flavor pairing of u and d quarks with 
additional single-flavor pairing of the remaining species ("2SC+1SC"), single-flavor 
pairing only ("ISC"), and three- flavor crystalline ("LOFF") pairing. 

These phases arise out of the competition between three factors: the strange 
quark mass, which increases the energy cost of s quarks, depressing their Fermi 
momentum; the requirement of electrical neutrality, which increases the number of 
u quarks to compensate for the smaller number of s quarks; the energy gain that 
quarks obtain by pairing. 

1. Single-flavor pairing (ISC). 

If the strange quark mass is large enough then charge neutrality requirements 
will force all three flavors to have Fermi momenta too far apart to allow any 
BCS pairing of different flavors. In this case, each flavor will have to self-pair 
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Figure 7: Pictorial representation of single flavor pairing ("ISC") in neutral quark 
matter. The requirement of electric neutrality and a nonzero strange quark mass 
forces the Fermi momenta of the three flavors apart. The red and green colors of 
each flavor pair in a color-antisymmetric channel. The blue s quark self-pairs in a 
color and flavor symmetric channel. We have not found an attractive channel in 
which the blue u or blue d could pair. Competing pairing patterns for the same 
arrangement of Fermi momenta are color-spin locking and the CFL-LOFF phase 
(see text). 

(Fig. 7) 1 . Two colors (red and green, say) can pair via gluon interactions 
in the (3a,3s,1,+)(C73) channel, and the blue s quarks can then pair in the 
rotationally invariant (65,35,0,+) (C7075). Note that these pairing patterns, 
being flavor symmetric, only exist for the single-gluon interaction, not the 
instanton interaction, so their strength is model-dependent. This leaves the 
blue u quarks and d quarks. We have not found a single- color and single-flavor 
channel that can pair massless quarks, so they will have to find an attractive 
channel outside the set that we have explored here. An alternative possibility 
is that for each flavor there is color-spin locked (CSL) condensation of all three 
colors. 

2. Two-flavor u-d and additional single flavor pairing (2SC+1SC). 

It is just possible that there is a range of m s in which 2SC u-d pairing can 
occur. m s must be high enough so that CFL is disfavored relative to 2SC, but 
low enough so that the u and d Fermi momenta are close enough to allow them 
to pair with each other. The energy accounting is complicated by contributions 
from the chiral condensate, which competes more or less strongly with different 

1 For a strange quark mass that is large (probably unphysically so), there will be a region in 
the quark matter part of the m s -fi plane where [i < m s , so there are no s quarks and the u and 
d Fermi momenta differ by some appreciable fraction of fi. This region is described by Fig 7 or 
Fig 8 but with the strange quark Fermi momentum set to zero, so no pairing of strange quarks. 
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Figure 8: Pictorial representation of 2SC+S with additional single-flavor pairing in 
neutral quark matter. This will only occur if the condensation energy of the 2SC 
pairing is strong enough to offset the cost of dragging the red and green u and d 
Fermi momenta away from the values dictated by electric neutrality and the strange 
quark mass, to a common value. 

quark pairing patterns. Current indications are that 2SC either does not occur, 
or survives only in a very narrow window [18, 19]. The pairings involved in 
2SC+1SC are depicted in Fig. 8. The red and green u and d quarks undergo 
2SC pairing. The blue quarks and the strange quarks follow the same pairing 
patterns as in the single-flavor pairing (ISC) case. Note that enforcing color 
neutrality will create a small 0(A 2 / fi) splitting between the Fermi momentum 
of the blue s and that of the the red and green s. This may prevent color- 
spin-locking for the s quarks. 

3. Three-flavor crystalline pairing (LOFF). 

An important competitor to the channels discussed here is the LOFF phase, 
in which flavors that cannot undergo conventional BCS pairing because their 
Fermi momenta are too far apart can nevertheless pair over part of their 
Fermi surfaces by forming pairs with net momentum. This breaks translational 
invariance, and leads to a crystal structure [24, 25]. One can imagine a type 
of CFL-LOFF pairing, in which the red and green u and d quarks, the red and 
blue u and s quarks, and the blue and green d and s quarks, all form LOFF 
pairs. A recent Ginzburg-Landau study [26] indicates that the condensation 
energy of the crystalline phase may be large enough to beat any single-flavor 
channel, although this result followed from extrapolating the Ginzburg-Landau 
analysis beyond its realm of validity. 

As well as evaluating the attractiveness of a large set of channels, we have 
calculated the angular dependence over the Fermi surface of the energies of the 
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quasiquark excitations associated with each of the j = 1 quark pair condensates. 
We see that the j — 1, rrij — ±1 condensates always have gapless quasiquarks at the 
poles of the Fermi sphere. The j = l,rrij = condensates have gapless regions for 
massless quarks, at the poles {ipC^ip) or around the equator (ipCao^ip), but if the 
quarks are massive then the the quasiquarks have a minimum gap of order m A / fi 
(Figs. 5,6). The gapless points or lines on the Fermi surface are likely to lead to 
interesting phenomenology (see below). 

Finally, we have surveyed the isotropic mesonic condensates, and find that 
the conventional chiral condensate is very heavily favored. The color-flavor-locked 
mesonic channel studied in Ref . [29] is strongly disfavored, as was noted in Ref . [36] . 
It has been argued [37] that taking into account the effect of the color-flavor-locked 
channel on the integral over instanton sizes makes the condensation energy come 
out attractive after all. 

7 The road ahead 

This work can be extended in many directions, and points to interesting questions 
in the theory of dense matter and compact star phenomenology. 
Dense matter: 

(1) In order to allow comparisons between the phases we studied and other 
competing possibilities such as the LOFF phase, it would be useful to calculate 
the free energy of the various ISC pairing patterns. (2) It would be valuable to 
extend our survey to include an even wider set of possible quark pair condensation 
patterns, such as color-spin locking. We would also like to find a channel in 
which light single-color single-flavor pairing can occur. It might be necessary to 
explore condensates where orbital angular momentum is introduced explicitly via 
derivatives: ijjCd^ip etc. These momentum-dependent condensates will not arise 
from the mean-field treatment of the simplified point-like interactions that we used. 
One would have to give the four-fermion vertex a non-trivial form factor, or go 
beyond the mean-field approximation. (3) A related area that is already being 
investigated is single color/flavor channels in the ultra- high-density region with 
the true QCD interaction vertex, including the hard-dense-loop-resummed gluon 
propagator [27, 15, 32]. In this regime one will also find momentum-dependent gaps. 
(4) We studied the attractive channels independently of each other, and independent 
of chiral condensation. It is clearly necessary to solve their coupled gap equations, 
and see how they compete/coexist. (5) We only worked at zero temperature. It is 
straightforward to generalize the gap equations to arbitrary temperature, and recent 
work indicates that the usual BCS relationship between the critical temperature and 
the gap parameter may be modified for these exotic condensates [28, 32]. 
Compact star phenomenology: 

Although many of the channels we studied have very small gaps and therefore very 
small critical temperatures, they will be phenomenologically relevant if they are the 
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best pairing option available for some of the quarks. Since the temperature of a 
compact star falls to tens of eV when its age reaches about a million years, pairing 
can suddenly occur in such channels late in the star's life, and the corresponding 
quasiquark excitations will suddenly become too heavy to participate in transport 
processes. We must remember, however, that some of the exotic channels are gapless 
at special 3-momenta. In those small proportion of their quasiparticles 

may continue to play a role in transport properties, even when the temperature is 
much less than the gap parameter. Some specific topics for further investigation 
are: (1) Develop the transport theory of j ^ condensates: they have gapless 
modes at points or along lines on their Fermi surface. What effect does this have 
on neutrino emission/absorption via URCA processes or otherwise, their specific 
heat, viscosities, conductivities, etc? A natural first step would be to write down 
an effective theory, which would contain the lowest quasiquark modes, Goldstone 
bosons arising from the breaking of rotational symmetry (which could be called 
"spin waves" by analogy with 3 He), and density waves. (2) The j = I condensates 
can carry angular momentum simply by aligning themselves in large domains, 
without involving any superfluid vortices, but it seems they will typically occur 
in conjunction with other phases that are superfluid. It would be interesting to see 
how the angular momentum is carried in this situation. (3) The j = 1 condensates 
are ferromagnets. It would be interesting to see if they could generate large magnetic 
fields such as those suggested for magnetars (10 14-15 Gauss). (4) The various j = 1 
condensates show completely different variation of the energy gap over the Fermi 
surface (Section 4). It would be useful to know how they behave when exposed 
to the nonisotropic external influences that are common in compact stars, such as 
magnetic fields or neutrino fluxes, and also in their coupling to external sources of 
torque, eg via electron-quasiquark scattering. 

Clearly there remain many interesting questions, both formal and phenomeno- 
logical, about the single color and/or single flavor color superconducting phases of 
dense quark matter. 
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A Calculational details 

To allow the possibility of quark pairing, we use 8-component Nambu-Gorkov 
spinors, 

-WW ™ 

with 

tt = {4>{p)^ T {-p)) . (A.2) 

In Minkowski space the inverse quark propagator for massive fermions takes the 
form, 

where 

A = 7o A t 7o . (A.4) 

The gap matrix A is a matrix in color, flavor and Dirac space, multiplied by a gap 
parameter also denoted as A, 

A£f° 6 = A(/i)C Q/ %T afe . (A.5) 

The relation between the proper self energy and the full propagator is, 

S " = V + E = (' - m + " 7 " „ + m °_ „ 7o)r ) + ( A ) <*<» 

where S*^ 1 is the inverse propagator in the absence of interactions. The gap is 
determined by solving a self-consistent Schwinger- Dyson equation for S. For a 4- 
fermion interaction modelling single gluon exchange, this takes the form 

£ = -6*G | *?LvtS(p)V+> (A.7) 

where V^ 4 is the interaction vertex in the Nambu-Gorkov basis. We study three 
interactions, the quark-gluon vertex 

A _( 1 ,X A /2 \ 

the quark-magnetic gluon vertex 

and the quark-instanton vertex, for which 

£? fc 7 = -GiG J -J^ (v&Sff(p)V*> + VlSff^V^)^ (A.10) 
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where 
and 



,2 



;i + 75) 







(i+ 75 ) t ; ' a ^ vr ~\ 



and V"r 



;i - 75) 







(A.11; 



(A.12) 



In the case of the i/jC^^i/j condensate for the full gluon interaction we obtain the 
gap equation, which after rotation to Euclidean space becomes 

dp d 3 P 4(A 2 + fi 2 +p 2 +p 2 ) 

(A.14) 
(A.15) 

(A.16) 



(2VT)* W 
where 

W = A 4 + ^ + (p 2 + p 2 ) 2 + 2A V + p 2 + p 2 ) - 2/i 2 (- Po 2 + P 2 ) 
The integral can be explicitly evaluated, 



2G f A , 
1 = — / dp 



7T Z 



v/A2 + (p + / x)2 y/A' + ip-rf 
The momentum integral can be performed analytically, giving 

r 2 



A = 2VA 2 -/i 2 exp ( A 2 ^ 2 3/i ] exp 



7T" 



for A < fj,. 



B Gap Equation Summary 



Here are the gap equations for the attractive channels. In the following, positive 
square roots are implied and we define p 2 = (p x ) 2 + (p y ) 2 - 



d\p\ 



d\p\, 



dp r dp z 



dp r 



B.l C75 and C gap equations 

G 



N- 



7T 



d\p\ 



\P\ 



+ 



dp z 



N 
N 
N 



4 

2 

3 
2 



(b.i; 



/ A 2 + (b|-/i) 2 v/A 2 + (H+/i) 5 
where N is a constant that differs for each interaction. 
Instanton 

Magnetic + Electric Gluon 
Magnetic Gluon 

The C channel produces an identical gap equation for both the full gluon and 
magnetic gluon interactions. The instanton interaction is not attractive in this 
channel. 
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B.2 C(To3 an d C<Jo375 gap equations 



G 

1 = N— I dp r dp 

IX 1 ' 



p r (S +p 2 r ) p r (£ -p. 



EE, 



with 



£ 2 = A 2 Pr 2 + /, 2 H 2 

A 2 + fi 2 + \p\ 2 ±2£ 



El 



where N is a constant that differs for each interaction. 



Instanton 
Magnetic Gluon 



N 
N 



(B.2) 



The Ccro375 channel produces an identical gap equation for magnetic gluon 
interaction. The instanton and the magnetic gluon interactions are not attractive 
in this channel. 



B.3 C(a i ±icr 02 ) gap equation 



N—f- I dp r dp z 



dp 



Pr jfJ? ~ (P°) 2 ~ pj ~ 2p°p z 

W 



Instanton 
Magnetic Gluon 



N = 2 
N= ± 



(B.3) 



where 

W = // + (~V) 2 + \p\ 2 ) 2 + 2A V - (p ) 2 -p\- 2p° Pz ) - 2^f + |p| 2 ) 
and N is a constant that differs for each interaction. 



B.4 C73 gap equation 

G 



N- 



7T 



dp r dp z 



Pr(£ +P 2 Z ) , Pr(S-p 2 z 



£E A 



with 



£E_ 



£ 2 = A 2 p 2 + /i 2 (H 2 + m 2 ) 

A 2 + fi 2 + m 2 + \p\ 2 ±2£ 



El 



where N is a constant that differs for each interaction. 

Magnetic + Electric Gluon N 
Magnetic Gluon N 



(B.4) 
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B.5 C7375 g a P equation 



This channel is not attractive for instantons and the magnetic gluon interaction is 
the same as for the massless C73 channel, i.e. 



1 = N 



G 



dp r dp z 



EE A 



EE. 



with 



E 2 = A 2 p 2 z + p 2 \p\ 2 

A 2 + fi 2 + \p\ 2 ±2E 



El 



(B.5) 



B.6 C(7i ±^72) gap equation 



N 



-iG 



7T" 



dp r dp z 



dp 



p r (2A 2 E^E 1+ + (m 2 + /i 2 - (p ) 2 + p 2 z )E 2 .E 2+ ) 

E^E l+ (AA A + 4A 2 (m 2 + /x 2 - (p ) 2 + p 2 z ) + E^E 1+ ) 

(B.6) 



with 



E\± = m 2 — (p ± p 



0\2 



Pz 

E 2 ± = m 2 - (fi±p ) 2 + \p\' 

where N is a constant that differs for each interaction. 

Magnetic + Electric Gluon 
Magnetic Gluon 



N 
N 



B.7 C7075 gap equation 
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1 = N- d\p\ 
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EE. 



with 



E 2 = A 2 H 2 + yU 2 (H 2 + m 2 
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A 2 + /i 2 



\p\ 2 ±2E 



where N is a constant that differs for each interaction. 

Magnetic + Electric Gluon N 
Magnetic Gluon N 



For m = this reduces to 



N 



2G A 3 
3^ ^A 2 + 



(B.7) 



(B.8) 
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C Orbital/spin content of the condensates 



In the non-relativistic limit it is meaningful to ask about the separate contributions 
of the orbital and spin angular momenta to the total angular momentum of the 
diquark condensates. We can find out by expanding the field operators out of which 
the condensates are built in terms of creation and annihilation operators, 



m = £ 



VE k 



u\k)a s kia e- ik * + v s {k)biJ k - 



(CI) 



Inserting the explicit momentum-dependent spinors in any basis allows the 
creation/annihilation operator expansions of the condesates to be calculated. 
In the Dirac basis, 



uf(k) = A 



uf (k) = A 



( 1 \ 



\B(h + ik 2 ) J 
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(C.2) 



.4 



E + m\ 2 
2m 



B 



E + m 



(C.3) 



Eq. (C.4) shows the result of performing such a calculation for the ipCip 
condensate, 



i 2 2 



pia —pj(3 ' pia —pjf3 



6 Sa 6 -»ifl)(Pl 



W2j 



pia —pjf3 ' u pia —pj[) 



S% (C4) 

where S^n is the color-flavor matrix which is symmetric under the interchange i ^ j 
(flavor) and a ^ (3 (color). A sum over momentum p should be performed on the 
right hand side. 

Once the operator expansions have been obtained it is a relatively simple 
procedure to obtain the angular momentum content by rearranging the terms and 
inserting the relevant spherical harmonics. It is important to include contributions 



+ 2 ( a lia a -pjp + V- P jf3)P3 
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from momenta k and —k together, since they involve the same creation/annihilation 
operators. For example, for the condensate in (C.4), 



p = k : — [aliad^pih - ik 2 ) - aL a -fe^(^i + ifa) + 2a lt^ 2 -kj pfa] (C.5) 
P=-k : - [-a^aljpih - ik 2 ) + a}_ kia a l kjp (k x + ik 2 ) - 2a 1 _ kia a 2 kjp k 3 ](C.6) 
Hia^-kjisih - ik 2 ) - al^a^jpih + ik 2 ) + 2a 2 kia a l _ kjf3 k 3 ] 

—k, i <-> j, a <-> /3 in the last line. The final result is 



E 
1 

1 



where we have relabelled k 
a sum over k, a, (3, i, j of 



E 



Upon inserting the relevant spherical harmonics and using standard arrow notation 
for the spins we obtain 



L7I ITT)rfl + 7! Iu)nl " 7!7! [ITi)+ liT)]Fl ° 



(C.8) 



which has precisely the correct Clebsch-Gordan structure to be interpreted as a state 
with orbital angular momentum 1 = 1, which gives an antisymmetric spatial wave 
function, and spin s = 1, which gives a symmetric spin wavefunction, combined to 
give j = 0. We write this as \l = 1a, s = Is). Applying this to all the condensates 
we studied, we can make a table of the particle-particle (as opposed to particle-hole) 
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content of each of them, 
4)C<r 03 A4> 
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These results are summarized in the "BCS-enhanced" column of Table 1. We can 
see explicitly that the ipC^o^^ip condensate has no particle-particle component in 
the massless limit, which is why it cannot occur at high density for the up and down 
quarks. This reflects basic physics: the condensate has spin zero, so the two spins 
must be oppositely aligned. But it is an "LR" condensate (see Table 1), so in the 
massless limit the two quarks, having opposite momentum and opposite helicity, 
have parallel spins. 



References 

[1] B. Barrois, Nucl. Phys. B129 (1977) 390. S. Frautschi, Proceedings of workshop 
on hadronic matter at extreme density, Erice 1978. 

[2] D. Bailin and A. Love, Phys. Rept. 107 (1984) 325, and references therein. 

[3] M. Iwasaki, T. Iwado, Phys. Lett. B350, 163 (1995); M. Iwasaki, Prog. Theor. 
Phys. Suppl. 120, 187 (1995) 



28 



[4] M. Alford, K. Rajagopal and F. Wilczek, Phys. Lett. B422, 247 (1998) [hep- 
ph/9711395]. 

[5] R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, Phys. Rev. Lett. 81, 53 
(1998) [hep-ph/9711396]. 

[6] M. G. Alford, Ann. Rev. Nucl. Part. Sci. 51 (2001) 131 [hep-ph/0 102047]. 
K. Rajagopal and F. Wilczek, hep-ph/0011333. T. Schafer and E. V. Shuryak, 
Lect. Notes Phys. 578 (2001) 203 [nucl-th/0010049]. D. K. Hong, Acta 
Phys. Polon. B 32 (2001) 1253 [hep-ph/0 10 1025]. S. D. Hsu, hep-ph/0003140. 
D. H. Rischke and R. D. Pisarski, nucl-th/0004016. 

[7] M. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. B537, 443 (1999) [hep- 
ph/9804403]. 

[8] J. Bardeen, L. Cooper, J. Schrieffer, Phys. Rev. 106, 162 (1957); Phys. Rev. 
108, 1175 (1957) 

[9] G. Carter, D. Diakonov, Phys. Rev. D60, 016004 (1999) 

[10] J. Berges and K. Rajagopal, Nucl. Phys. B538, 215 (1999) [hep-ph/9804233]. 

[11] T. M. Schwarz, S. P. Klevansky and G. Papp, Phys. Rev. C 60 (1999) 055205 
[nucl-th/9903048]. B. Vanderheyden and A. D. Jackson, Phys. Rev. D 62 
(2000) 094010 [hep-ph/0003150]. F. Gastineau, R. Nebauer and J. Aichelin, 
Phys. Rev. C 65, 045204 (2002). B. O. Kerbikov, hep-ph/0106324. M. Huang, 
P. f. Zhuang and W. q. Chao, Phys. Rev. D 65 (2002) 076012 [hep-ph/01 12124]. 
M. Oertel and M. Buballa, hep-ph/0202098, Nucl. Phys. A 703 (2002) 770 [hep- 
ph/0109095]. 

[12] S. Hands and D. N. Walters, hep-lat/0209140. 

[13] D. Son, Phys. Rev. D59, 094019 (1999). 

[14] T. Schafer, F. Wilczek, Phys. Rev. D60, 114033 (1999). R. Pisarski, D. Rischke, 
Phys. Rev. D61, 074017 (2000). D. Hong, Nucl. Phys. B582, 451 (2000) [hep- 
ph/9905523]; D. Hong, Phys. Lett. B473, 118 (2000) [hep-ph/9812510]. R. 
Pisarski, D. Rischke, Phys. Rev. D61, 051501 (2000). D. Hong, V. Miransky, 
I. Shovkovy, L. Wijewardhana: Phys. Rev. D61, 056001 (2000); erratum Phys. 
Rev. D62, 059903 (2000). S. Hsu, M. Schwetz, Nucl. Phys. B572, 211 (2000). 

[15] W. Brown, J. Liu, H. Ren: Phys. Rev. D 61, 114012 (2000) [hep-ph/9908248]. 
Phys. Rev. D 62, 054016 (2000) [hep-ph/9912409]. Phys. Rev. D 62, 054013 
(2000) [hep-ph/0003199]. 

[16] T. Schafer, Nucl. Phys. B575, 269 (2000) [hep-ph/9909574]. 



29 



[17] I. Shovkovy, L. Wijewardhana, Phys. Lett. B470, 189 (1999). 

[18] M. Alford and K. Rajagopal, JHEP 0206 (2002) 031 [hep-ph/0204001]. 

[19] A. W. Steiner, S. Reddy and M. Prakash, hep-ph/0205201. 

[20] T. Schafer and F. Wilczek, Phys. Rev. D 60 (1999) 074014 [hep-ph/9903503]. 

[21] M. Alford, eConf C010815 (2002) 137 [hep-ph/0110150]. 

[22] P. F. Bedaque and T. Schafer, Nucl. Phys. A 697 (2002) 802 [hep-ph/0105150]. 

[23] T. Schafer, Phys. Rev. D 65 (2002) 094033. 

[24] A. I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 47, 1136 (1964) [Sov. 
Phys. JETP 20, 762 (1965)]; P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 
(1964). 

[25] M. Alford, J. Bowers and K. Rajagopal, Phys. Rev. D 63, 074016 (2001) [hep- 
ph/0008208]. 

[26] J. A. Bowers and K. Rajagopal, Phys. Rev. D 66, 065002 (2002) [hep- 
ph/0204079]. 

[27] T. Schafer, Phys. Rev. D62, 094007 (2000) [hep-ph/0006034]. 
[28] M. Buballa, J. Hosek and M. Oertel, hep-ph/0204275. 

[29] C. Wetterich, Phys. Rev. D 64, 036003 (2001) [hep-ph/0008150]. J. Berges and 
C. Wetterich, Phys. Lett. B 512, 85 (2001) [hep-ph/0012311]. 

[30] B. Barrois, "Nonperturbative effects in dense quark matter", Cal Tech PhD 
thesis, UMI 79-04847-mc (1979). 

[31] N. K. Glendenning, Phys. Rev. D 46, 1274 (1992); N. K. Glendenning, Compact 
Stars (Springer- Verlag, 1997); F. Weber, J. Phys. G. Nucl. Part. Phys. 25, R195 
(1999). 

[32] A. Schmitt, Q. Wang and D. H. Rischke, nucl-th/0209050. 

[33] M. Alford, J. Berges and K. Rajagopal, Nucl. Phys. B558, 219 (1999) [hep- 
ph/9903502]. 

[34] T. D. Fugleberg, hep-ph/0206033. 

[35] R. D. Pisarski and D. H. Rischke, nucl-th/9907094. 

[36] J. Berges, Phys. Rev. D 64, 014010 (2001) [hep-ph/0012013]. 

[37] C. Wetterich, Phys. Lett. B 525, 277 (2002) [hep-ph/0011076]. 



30 



